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Abstract 

The hadronization process for quarks combining into two mesons, 
qq — ► MM' at temperature T is described within the SU(3) Nambu- 
Jona-Lasinio model with finite current quark masses. Invariant matrix 
elements, cross-sections and transition rates are calculated to leading 
order in a 1/N C expansion. Four independent classes, ud, us, uu and 
ss — * hadrons are analysed, and the yield is found to be dominated by 
pion production. Threshold behaviour is determined by the exothermic 
or endothermic nature of the processes constituting the hadronization 
class. A strong suppression of transition rates is found at the pio- 
nic Mott temperature Tm-k = 212 MeV, at which the pion becomes a 
resonant state. The mean time for hadronization is calculated to be 
2-4fm/c near the Mott temperature. The calculation of strangeness 
changing processes indicates that hadronization accounts for a 1% in- 
crease in the absolute value of the kaon to pion ratio at T = 150 MeV. 
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I. INTRODUCTION 



One of the current outstanding problems facing nuclear and particle theoreticians 
today lies in understanding the phase transition from a constituent quark and gluonic 
matter to that of observed hadronic matter. Since it is expected firstly that such a 
deconfined state of quarks and gluons exists, and secondly that it may be observed 
via Pb + Pb and other heavy ion collisions already and also to be undertaken at 
CERN, an adequate description of the hadronization process that leads one from 
the quark and gluon degrees of freedom inherent to QCD to the observed hadron 
spectrum is really required. The present state of the art (for an overview over the 
field see the proceedings of the Quark Matter conferences (l],[2j or the book ||) 
is characterized either by phenomenological approaches for the parton to hadron 
transition or computer codes based on string phenomenology ||. 

In this paper, we place our emphasis on describing the hadronization process 
of quarks and antiquarks into two mesons, within a microscopic, field theoretical, 
and non-perturbative framework, which is carried through at finite temperatures. 
In itself this is a demanding project that cannot yet be performed directly starting 
from the QCD Lagrangian, and thus the price that we have to pay is in making 
a choice of a model Lagrangian. To this end, we invoke the Nambu-Jona-Lasinio 
(NJL) model [7|-[| in its SUl(3) x SUr(3) version, which has been constructed to 
display the same internal symmetries as QCD itself. This model is known to provide 
a good description of the static properties pertaining to both the nonstrange and 
strange meson sector at zero temperature, and it allows for a transparent description 
at finite temperatures, explicitly displaying the chiral symmetry restoration phase 
transition at a critical temperature of about 200 MeV in the chiral limit. Although 
we are aware of the deficiencies of the model — lack of confinement as well as 
non-renormalizability in a strict field theoretic sense — it nevertheless is possible 
to construct a comprehensive physical picture from the calculations performed to 
date both at T = and at finite temperatures and densities. For definiteness, it is 
useful to argue from the chiral limit. Assuming that the chiral and deconfinement 
phase transitions occur at the same temperature — an assumption supported by 



lattice gauge simulations [10] — one can build a heuristic picture on the basis of the 



NJL model. At T > T c , hadronic states are unstable, and the system is dominated 
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by interacting quarks with constituent mass m — 0. For T < T c , where chiral 
symmetry is broken, m^O, m n = 0, m a = 2m, and a plasma of quarks and mesons 
is formed. While the appearance of quarks for temperatures T < T c is an artifact of 
the model due to the absence of confinement, the situation for T > T c is realistic, 
except for the lack of explicit gluonic degrees of freedom. As such, however, the 
process of hadronization may be reliable around the critical temperature, where the 
shortcomings of the model may not be too severe. These arguments can be extended 
to the physical case of non-zero current quark masses, which we in fact study here. 

The purpose of this paper is thus to investigate all hadronization processes of 
the type qq — > MM', where M represents a meson in the non-strange or strange 
sector, and here q = u,d, s. In particular, we evaluate scattering amplitudes and 
cross-sections, and study the hadronization rates for all these processes, which in 
turn are necessary elements for constructing a dynamical non-equilibrium transport 
theory for this model Lagrangian [[□]]. This should aid us in understanding the role 
that chiral symmetry plays in dynamical processes. Our approach follows in part, 
and develops further the calculation performed in SU{2) in the chiral limit described 
in the letter jg. 

In our study of the SU (3) sector, we classify hadronization processes according 
to the incoming quark and antiquark, finding a total of four independent classes 
under the additional assumption of SU{2) isospin symmetry, where m u = m^. We 
thus consider the processes with ud, us, uu and ss as incoming pairs, which we also 
list. The mesons M and M' considered in the process are tt, K and m which are the 
stable ones with respect to the strong interaction. Feynman diagrams for the scat- 
tering amplitudes pertaining to these reactions have a generic form which describes 
all qq —>■ MM' . We select the diagrams according to an expansion in the inverse 
number of colors in the model, 1/N C , that is commensurate with the Hartree approx- 
imation for the self-energy and the random phase approximation for the scattering 
amplitude |T3HT5|] . The variations in applying each of the Feynman diagrams to the 
different scattering channels arise from differing exchanged mesons, quark masses 
and flavor factors. We choose to illustrate our calculation by analytically construct- 
ing the formulae for the cross sections of the processes ud — > 7r + 7r°, us — > tt + K° and 
uu — > 7r + 7r~. The first of these explicitly includes a u channel exchange, while the 
second has unequal masses for all incoming and outgoing particles. The final process 
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illustrates the role of the mixing of the scalar mesons in the £77(3) Lagrangian due 
to the t'Hooft term and the differing up, down and strange quark masses. 

All quantities are calculated as a function of temperature and baryonic chemical 
potential. It turns out that a cardinal role is played by the Mott temperatures, Tm^, 
Tm<ji Tmk, Tm v , that are defined to be the temperatures at which the masses of the 
respective mesons tt, a, K, rj, are equal to the sum of the masses of their constituents. 
Note that these four mesons are bound states at T = 0. At the Mott temperatures, 
the respective meson-quark couplings go to zero, and this in turn influences the 
transition amplitudes and transition rates in that they are suppressed. In particular 
the pionic Mott temperature plays the dominant role, driving the total transition 
rates (almost) to zero at this point. We comment that the use of a finite temperature 
explicitly breaks Lorentz invariance in our calculation. We do not attempt to resolve 
this problem here. Quantities are simply calculated with respect to the rest frame 
of the medium. 

In order to make a connection with transport theory, we calculate transition 
rates, which are constructed by multiplying the total cross section of a particular 
process with the relative velocity. Our numerical results indicate that in a per- 
channel calculation, the t and u channel exchange of any one process dominate 
over the s channel exchange, except perhaps in an energy range where a resonance 
is present. The t and u channel graphs, if present in a process, give cancelling 
contributions at threshold when the flavor factors are equal. In dealing with the 
four independent hadronization classes, we are directly able to identify the leading 
contributions at all energies. These, with the exception of the ss hadronization class, 
are dominated by pion production. The ss processes on the other hand are found 
to be dominated by kaon production. In all these calculations we can confirm that 
the exothermic reactions are always dominant over the endothermic ones, as far as 
they are present. 

We also examine the energy averaged transition rates. These quantities, mul- 
tiplied by the density of incoming quarks, can be interpreted as the inverse of the 
hadronization time. We find this to be of the order of r = 2-4fm/c in the tem- 
perature range 160-200 MeV, rising rapidly as one moves towards the pion Mott 
temperature Tm-k, or to lower temperatures. This indicates that hadronization oc- 
curs preferentially over this given range of temperatures within our model. 
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Finally we have also examined processes that change the number of s plus s 
quarks, like e. g. ud — > K + K°. We find that the hadronization process generates 
an enhancement of 1% in the absolute value of the ratio of kaons to pions at T — 
150 MeV. This is to be compared with the experimental value of (15.4 ± 0.8)% for 
the Kj it ratio found for S + S collisions, which itself constitutes an enhancement of 



9% over the N + N value of (6.5 ± 1.1)% p 



In order to facilitate the technical evaluation, we decompose all quantities into 
fundamental integrals containing one, two or three denominators, and which arise for 
instance from the self-energy, polarization and the three meson vertex. In principle, 
these functions can be calculated for arbitrary values of the differing quark masses 
and associated differing chemical potentials. In this way, general expressions can 
be constructed if one utilizes the modular forms. For actual calculational purposes, 
however, we will set all chemical potentials to zero, and study temperature and 
energy dependencies here, giving some relevant analytical results for these integrals 
in an appendix. 

This paper is structured as follows: In Section ||, we discuss the general SUl(3) x 
SUr(3) Lagrangian that we use, and give the relevant functions associated with the 
mass spectrum evaluation that we require. In particular the rj-rf mixing is detailed, 
and the scalar resonance sector, which plays the role of intermediate states in the 



hadronization. In Section |TTT| , the possible hadronization processes are classified, 
and the explicit examples ud — > 7r + 7r°, us — > tt + K° and uu — > 7r + 7r~ are discussed. 
Numerical results for the transition rates for each hadronization class are given in 
Section [IV]. Strangeness changing processes are also analyzed here. We summarize 
and conclude in Section [Vj. 

II. PROPERTIES OF MESONS IN SU(3) 

A. General Considerations 

This section serves to introduce our Lagrange density and notation, commencing 
with the three flavor Lagrangian 

3 

£ = 53 - mofWf 
f=o 
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a=0 



-K [det V(l + 75)^ + det ip(l - 75)^ 



(1) 



where G and are dimensionful coupling strengths. The term containing G dis- 
plays U(3) x U (3) symmetry, while the determinantal term controlled by K breaks 
this down to SU{3) x SU(3). Flavor and color indices have been suppressed for 
convenience in the interaction terms. However, the flavor indices have been explic- 
itly included in the first term, since the current quark masses m /, that themselves 
explicitly break the SU(3) x SU(3) symmetry, are regarded as distinct. For gen- 
eral reviews on the three flavor version of the NJL model, the reader is referred to 



It is useful to convert the determinantal term in Eq. (|1|) into an effective two- 
body term in the mean field approximation. This follows on contracting out one pair 
of quark and antiquark fields and dividing this result by two. One may recombine 
the result with the existing two body term that is controlled by G to find the effective 
Lagrangian to be 



+ir 30 (^A 3 ^)(^AV) + ^(^a^x^aV) 
+ir 03 (^AV)(^A 3 ^) + K+^AVX^A 3 ^) 

+K 8 - (^AV)(^AV) + ir+(^7 5 AVX^75AV) 
+K - 8 (^A ( VX<M 8 VO + K+(^75AVX^75AV) 
+^ 83 (^AVX<M 3 ^) + ir+^AVX^A 3 ^) 
+K 38 (^A 3 ^)(^AV) + ^3 + 8(^75A 3 ^)(^7 5 AV) (2) 



Refs. 
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+ J2 [K-(^Y + iC(^7 5 A» 2 



a=0 



with the effective coupling constants 




(3a) 




(3b) 
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Kf = Kf = G± ^N c Kitr 7 S d {x, x) (3c) 



Kq =Kf = G± -N c Kitr 7 S u (x, x) (3d) 



Kf = G ± ^N c K{2itr^S u {x, x) + 2iti ry S d (x, x) - tr 7 S s (x, x)) (3e) 
= = T^iV c #(ztr 7 S u (x, x) - iti^x, x)) (3f) 



/2 

K± = K± = ±y^N c K{itr 7 S u (x, x) + itr 7 S d (x, x) - 2itr 7 S s (x, x)) (3g) 



tf± = Kg = ±^-j=N c K(itr^S u (x, x) - itr^S d (x, x)) . (3h) 

In Eq. (|3|), tr 7 refers to the spinor trace alone, while is the diagonal quark 
propagator for a given flavor /, which, in the imaginary time formalism for finite 
temperatures can be written as 

P7 J {2ny j (iuj n + n f ) -jp-rrif 

Here the Matsubara frequencies are fermionic, iuo n = (2n + 1)tt/(3, with n = 

0, ±1, ±2, ±3, . . ., and fif is the chemical potential for a quark of flavor /. 

As can be seen from Eqs. (0, [|), the nondiagonal coupling constants Kq 3 , Kq 8 
and K% 8 give rise to the Tr°-r]-r]' mixing. If one assumes SU(2) isospin symmetry, 

1. e. m u = ma, then Kq 3 and K 38 vanish identically, with the consequence that the 
7T° decouples from the r? and rf . We make this assumption in what follows, using 
the generic label q for u, d, and always writing s explicitly. 

Writing 



itr 1 S f (x,x) = -—A(m f ,/j, f ) , (5) 



where 



M = }f^j^ ^ (6) 

\p\<A 



7 



(with Ej = p 2 + rrij) denotes the first loop integral (cf. Appendix one can easily 
derive the coupled gap-equations 



GN C A , , KN 2 C 

m u = m 0u —m u A(m u , fx u ) + -—-m d m s A(m d , fx d )A(m s , fi s ) (7a) 

7T Z 87T 4 



GN c A , x KNl 

m d = m od —m d A(m d , /i d ) + ——m s m u A(m s , n s )A(m u , fj, u ) (7b) 

1X Z 07T 4 



m , = mo , _ ^ m ,A(m„ „) + f? m „ m ^(»„, ^A (md , , d) (7c) 

from the mean field or Hartree approximation to the self-energy, and which deter- 
mine the physical quark masses. 

B. Pions and Kaons 

Imposing the degeneracy condition m u = m d , the determination of both the pion 
and kaon masses follows similarly to the standard approach taken in the two flavor 
model M. This comes about, since the off-diagonal coupling strengths Kq 3 and K^ 8 
are identically zero, and the 7r° is decoupled from the rj and rj . In this limit, the 
tx ± 1 7T° become degenerate. Concomitantly, Kf = K^, which means physically that 
neutral and charged kaons have the same mass. 

The quark-quark scattering amplitude is calculated in the random phase ap- 
proximation (RPA, see Fig. |I|), which yields the result 

MJk , k) = (8) 

1 - 4^+n^ , fc) 1 ; 

for the pion scattering amplitude and 

M K (k , k) = — * ^ (9) 

1-AKiXl^k) 

for the kaon one. Here Il p (/co, k) represents the irreducible pseudoscalar polarization 
(see Fig. 0), which, in the finite temperature Matsubara formalism, depends on ko 
and k separately, because the medium breaks Lorentz invariance. For arbitrary 
flavors, the irreducible polarization is given by the analytic continuation of the 
imaginary time form 
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-in 



1„ r d 3 p (iuj n + fj, 1 )(iu n -iu m + fi 2 ) ~p(p-k) -171^2 



d 3 p 



tr- 



iS fl (iu n ,p)ij 5 iS f2 (iu n - iv m ,p- k)i~f 5 



(10) 



where Ei = y p 2 + m\ , E 2 = y (p — k) 2 — m|, and the mesonic Matsubara frequen- 
cies iv m = 2m7r//3 are even, m = 0, ±1, ±2, . . .. Note that the irreducible polar- 
ization in this definition does not contain any flavor factors from the Gell-Mann 
A-matrices. These factors are incorporated explicitly in Eqs. (§, ||) as multiplicative 
cofactors of the coupling and polarization. 

We may decompose U p (ko, k) in terms of the function A(m, fi) already defined 
in Eq. (ffi) as 



n p (fc ,fc) 



N 



A[m^nx) + A(m 2 ,fi 2 ) 
+ ((mi - m 2 f - (k + fa - fi 2 ) 2 + A; 2 ) B (k, m 1 ,fi 1 ,m 2 , /i 2 , k ) 

introducing the second loop integral (cf. Appendix [A]) 



B (k,m 1 ,fi 1 ,m 2 ,fi 2 ,iu m ) 



(12) 



16tt 2 

P 



d 3 p 



1 



Ip1<a 



(27r) 3 ((iu) n + fii) 2 - Ef) ((iu n - iu m + n 2 f 



El 



(Note that in Eq. fllTD, an analytic continuation has been performed.) We comment 
that, due to rotational invariance, Bq does not fully depend on k, but only on \k\. 

In terms of the modular integrals, A and Bo, polarizations required can then be 
explicitly given for arbitrary chemical potential and temperature as 

Nr. 



n£(fc , k) 



nf ff (fc ,fc) 



N 



8tt 2 



2A(m q , fig) - k B (k, m q , fi q , m q , fi q , k ) 



2A(m s , fig) - k 2 B {k, m s , /i s , m s , /i s , k ) 



n£(*b, k) 



N c 


A( 


8tt 2 


+ [(m 


q — 



(13) 
(14) 
(15) 



The pion and kaon masses are determined according to the dispersion relations 
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(16) 



;i7) 



while effective couplings can be identified from the pole approximation forms 



M v (ko, k) 



-q 2 - 



k 2 — m 2 



to be 



M K (k , k) 



~9Kqs 



k 2 — m 2 K 



9-K. 



2 i anf ? (fc ,o) 



Q<1 



2m„ 



dkn 



k =m n 



(19) 



(20) 



9k 



1 0n£(fc o ,O) 



2m 



K 



dk 



k =m K 



(21) 



At low temperatures, Eqs. ( [Tq , 17) have bound state solutions with < 2m g 
and m^ < m q + m s . In this case, Eqs. (16| , |17"D are real equations. At higher 
temperatures, the polarization functions become complex functions with complex 
solutions for the meson masses, that we may write as p 



(22) 



m K ^m K - -Y K . (23) 

Denoting the Mott transition temperatures Tm-k and Tmk as the temperatures at 
which m n = 2m q and mx = rn q + m s respectively, one thus has the physical circum- 
stance that at temperatures larger than T Mn or T MK , respectively pions or kaons 
become resonances with finite widths due to the available decay channels into two 
quarks. From Eqs. (|20| , |2l|) , it follows that the quark-meson couplings also become 
complex in this case. 
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C. rj and rf 



Because of the mixing terms occurring in Eq. (§), the calculation of the rj and 
rf masses and couplings is somewhat more involved. The scattering amplitude is 
nondiagonal in this sector with entries Moo, Mg$, and Mq§ = M$o. Within the RPA, 
that is still expressed diagramatically in Fig. |], it can be calculated in matrix form 
to be § 

M = 2K + (1-2U P K + )- 1 , (24) 
where now K + and Il p are the 2x2 matrices 

,\ I K ' ,r K » H I (2r,) 



The Kf on the right hand side of Eq. (^5|) are as defined in Eqs. (|3|), while the 
pseudoscalar polarization functions are the linear combinations 

n P = \ (2l£ + nf f ) (27a) 



2 



n 8 p = 3 k + 2n y ( 27b ) 

< = <= 2 -^(K-^ s ) (27c) 



of the functions IT^ and IT^, that can be evaluated via Eqs. ([13], |14|). In Eqs. (^4| 



to (p7|), we have dropped the argument (k 0} k) for convenience. A determination of 
the masses and coupling strengths can be made on forming the inverse of the matrix 
M. We abbreviate this as 

M" 1 = — ^— r ~~ i (2? 

2detK+ 1 " 1 

with 
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A = K+ - 2IF°det J FsT + (29a) 

- B = K+ + 2U^ 8 detK + (29b) 

C = K+ - 2Tl£detK + . (29c) 
Following P|, we introduce the diagonal forms M v and M v ' via 

where it is a simple matter to verify that 



M' 1 = A + C - y (A - C) 2 + AB 2 (31a) 



M-, 1 = A + C + yf(A - C) 2 + AB 2 (31b) 



c 2 + s 2 = 1 (32a) 



c 2 -s 2 = C A (32b) 
'(A - Cf + AB 2 



-2B 

2cs = . (32c) 

y/(A - C) 2 + AB 2 

The masses of the rj and rj' mesons can now be determined via the condition 

M- 1 (m rn 6) = (33) 

Mj 1 (m v ,,6)=0 . (34) 

For meson masses that lie below the quark-ant iquark mass threshold, these equa- 
tions are again real. In practice, actual numerical determinations place the rj mass 
below this threshold at T = 0, and at a specific Mott transition temperature Tm^, 
this becomes a resonance, in much the same way as occurs for the pions and kaons. 
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The rf meson, however, is distinguished by the fact that it always lies above the 
quark-antiquark threshold, and it therefore is a resonant state at all temperatures. 
Since this is an artifact due to the lack of confinement, we regard this feature 
with some scepticism. We nevertheless calculate its contribution to the relevant 
hadronization cross sections and transition rates, and find that it is negligible, so 
that it can be safely discarded. 

To calculate the coupling constants from Eq. (p8|) , we express M directly as 



M 




(35) 



with 



V = (AC -B 2 )/ det K + 

= 1 - 2tr(rLftT + ) + 4 det n det K + . 



(36) 



At (k , k) = (m v , 0), V — 0, and in the usual fashion, we make a pole approximation 
giving the form 



M 



From this, one obtains 



4m r 



1 



3D 
dko 



ko=m v 



C -B 



k 2 -ml\-B A 



(37) 



with 



M. 



ab 



9ar]9brj 



2 

9o v 



2 

98r/ 



90ri98ri 



k 2 


— m 2 


Arrive 


dD 
dko 


ko=m v 




dP 
dko 


ko=m v 


4m v B 



dP 




dk 


k =rn v 



(38) 



(39) 



(40) 



(41) 



From these coupling constants we may calculate 
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9m = y ^9o n + ^9s v (42) 



2 2 

9 V ss = y g^o?? - y^sr, • (43) 

The coupling strengths g v > qq and g v ' sS may be evaluated in the same fashion. The set 
of couplings g vqq , g vsS and g v 'qq, g v 'ss enter directly into the cross section calculations. 



D. Scalar Resonances 

Since scalar mesons may occur as possible intermediate resonance structures in 
the hadronization cross sections, we are also interested in obtaining their masses. In 
contrast to the two flavor model, in which only the a is present, we deal here with 
nine scalar resonances: three ov's, which are the scalar partners of the pions, four 
or-'s, being the scalar partners of the kaons, and the a and a', that are associated 
similarly with the r) and rf . As occurs in the pseudoscalar case, for m u ^ ^ m s 
we have mixing between the a, a' and the neutral a n , the latter decoupling from 
the former if SU(2) isospin degeneracy, m u = ma, is imposed. 

The same techniques applied in Sees. II B| and [II Q can now be directly applied 



to the scalar resonances, with two changes: (i) one has to replace the coupling 
constants Kf of the previous section by K~ given in Eq. (|3|) and (ii) the pseudoscalar 
polarizations are replaced by their scalar counterparts: 

Rqq(ko,k) = -^ 2A(m g ,/j, q ) + (Aml~k 2 )B (k,m g ,ij, g ,rn g ,jjg,k ) (44) 



U^(k , k) = [2A(m s , /i s ) + (Am 2 - k 2 )B (k, m s , /i s , m s , fi s , fc )J (45) 



Rgsiko, k) = -— c - A(m q , fi 9 ) + A(m s , /x s ) (46) 



+ {{mq + m s ) 2 - (k + /J g - /j, s ) 2 + A" 2 ) B (k, m q , fi q , m s , fi s , k 

which are derived from the same graph as in Fig. 0, dropping the 275 factors at the 
vertices. 
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E. Numerical Results 



For our numerical calculations, we employ the parameter set mo q = 5.5 MeV, 
m 0s = 140.7 MeV, GA 2 = 1.835, KA 5 = 12.36 and A = 602.3 MeV, that has been 
determined on fixing the conditions m n = 135.0 MeV, tuk = 497.7 MeV, = 
957.8 MeV and f n = 92.4 MeV, while mo q is fixed at 5.5 MeV. The reason why we 
fitted the mass of the rj instead of the rj has a purely technical origin. This parameter 
set gives an 77 mass of m v = 514.8 MeV, which compares reasonably well with the 
physical value = 548.8 MeV. Although we have developed the general formalism 
to include the case of finite chemical potentials, we confine ourselves to [L q = fi s = 
in what follows. 

Figure ^| shows the temperature dependence of the constituent quark masses. At 
T = 0, we find m q = 367.7 MeV and m s = 549.5 MeV. At temperatures around 
200 MeV the mass of the light quarks drops to the current quark mass, indicating a 
washed out crossover from the chirally broken to approximately chirally symmetric 
phase. The strange quark mass starts to decrease significantly in this temperature 
range, but even at T = 300 MeV it is still a factor of two away from the strange 
current quark mass. 

Figure |] shows the temperature dependence of the pseudoscalar meson masses 
at n = 0. For comparison, the curves of 2m q and m q + m s are also indicated. At low 
temperatures, the meson masses are approximately constant. The crossing of the n 
and rj lines with the 2m q line indicates the respective Mott transition temperature 
for these particles, Tmk and Tm v - One observes that Tm v < 7W, the absolute values 
are Tm v = 180 MeV and Tmk = 212 MeV. For temperatures higher than Tm-k-, Tm v , 
respectively, the ir and rj become resonances and their masses increase. Similarly a 
Mott transition temperature Tmk for the kaon modes can be defined at the point 
where mx meets m s + m q . This is also indicated in the figure. One can also see 
from this plot that Tm-k and Tmk are approximately equal, with Tmk = 210 MeV. 

Figure [5] shows the absolute values of the pion and kaon coupling constants. 
A striking behaviour is observed at the Mott temperature. There the polarization 
displays a kink singularity, which can also be seen in the meson masses. Technically 
this results in the coupling strengths approaching zero for T — > Tm from below. This 
behaviour differs markedly from the behaviour of the couplings when evaluated in 
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the chiral limit. In that case, the coupling strength is always different from zero 
as one approaches the transition temperature from below. Note that this strong 
deviation from the chiral limit behaviour may have extreme consequences for results 
that depend strongly on this function, such as the cross sections. Since we regard the 
physical situation as being non-chiral, we investigate this situation only, and draw 
conclusions accordingly. The 77-quark couplings have not been explicitly shown; 
they display a qualitatively similar behaviour. 

The numerical calculations of the scalar mass spectrum are shown in Fig. |6]. 
At T = 0, we find m a = 728.9 MeV, m a , = 1198.3 MeV, = 880.2 MeV and 
m aK = 1050.5 MeV. For comparison, the double constituent quark mass is also 
shown in the figure. All these mesons are unstable over the entire temperature range, 
except for the a, for which the mixing results in a mass slightly below 2m q at T = 0. 
Although the difference is small compared to the standard two flavor model, it has 
the qualitative effect of making the a a stable particle for temperatures up to its 
Mott temperature Tmu = 165 MeV. As in the two flavor model, we obtain m a ~ m n 
above the pion Mott temperature, as is expected from symmetry requirements. 

The results presented in this section compare well with those of Ref. [|9] . 



Since a large number of hadronization processes are available to the light and 
strange sector quarks, it is useful to introduce a classification scheme to simplify 
the task of bookkeeping. The cross sections are classified according to the incoming 



fn u = nidi we have isospin symmetry and charge conjugation, leading to the relations 



III. HADRONIZATION CROSS SECTIONS 



A. Classification of Hadronization Processes 



quarks and include all exit channels Since we work in the approximation 



°ud — a du 




®us @ds O "su O "sd 



(47b) 
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using an obvious notation. Together with a sS , we thus have four independent classes 
of hadronization cross sections, cr u j, a uS , a ua and a sS . For these four classes, we 
determine the hadronization processes that are not forbidden by charge, strangeness 
or isospin conservation. These conservation laws lead to the processes that are listed 
in Table |. 

For the outgoing channels in Table |, further symmetry relations hold: 

&us-,tt+k° = 2&us-+ir K+ (48a) 

o"ud-»7r+7) = 2ow^7r°r) (48b) 
Vud-*K+rf — 2o u ^ 7r o 7? / (48c) 
Css-^K+K- = a ss^>K°K s (48d) 

Pss-tTT+TT- = 2<r s j_ >7r o 7r o , (48e) 

as a consequence of flavor algebra. In total, there are nineteen independent cross 
sections that we calculate. 



B. Feynman Graphs 

The Feynman graphs that we consider, have the generic forms shown in fig. |7]. 
This choice of diagrams is in keeping with the evaluation of the Hartree diagram 
to determine the gap equation and the random phase approximation for the po- 
larization. Together this selection constitutes a consistent expansion in the inverse 



number of colors, 1/N C |I3|-|T5[|. We stress that this is not an expansion in the 
coupling strength, but is rather a non-perturbative expansion. Note that the calcu- 
lation of the transition amplitudes is complicated with respect to the SU(2) isospin 
symmetric case by the fact that here each fermion line carries a flavor dependent 
mass, while the meson lines in turn carry differing masses also. We therefore at- 
tempt to retain as general a formalism as possible in what follows for the transition 
amplitudes and we specify the parameters later. 
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1. s channel 



The s channel exchange diagrams have the form 

- iM s = v(p 2 )u(p 1 )5 ClC J s iV(p 1 + p 2 )T(p 1 + p 2 ; p^)ig\ig% , (49) 

where g\ , g 2 are meson-quark coupling strengths for the outgoing mesons and f s is a 
flavor factor. The momenta of the incoming particles are pi and p 2 . For the outgoing 
particles, we assign momenta p 3 and p^. T> stands for the scattering amplitude of 
the virtual scalar meson. It can either take the form appropriate for a scalar meson 
that corresponds to Eqs. (|| [|), if the incoming quarks have different flavor, or it can 
be constructed as the sum over mixing terms according to Eq. ([24]), in all cases with 
Kf replaced by K~ , and IT P by II . The symbol F describes the three meson vertex 
contribution to the diagram. Its general form in the imaginary time formalism is 
(cf. Fig. |) 

-> i r d^q 

T(w m , k; ia h p) = ~N C - ^ J -^-^ (50) 



xtr„ 



iS h (iu n , q)ij 5 iS h (iuj n - ia h q- p)ij 5 iS h (iuj n - iu m , q- k) 



where (iu m , k) is the four momentum of the incoming scalar meson and (iai,p) the 
four momentum of one of the outgoing mesons. It is once again understood that 
the complex meson frequencies are to be analytically continued at the end of the 
calculation. Taking the spinor trace in Eq. ( pUD leads to the form 

T(iu m ,k;iai,p) = J ( 51 ) 

A — B 

x 



with 



[{iu n + /ii) 2 - El] [{iu n - iai + fi 2 ) 2 - El] [{iu n - iv m + /i 3 ) 2 - Ef\ 



A = m^q{q — p) — m 2 q(q — k) + rri\{q — p){q — k) + m\m 2 m^ (52) 

B = m 3 (iu n + ^i){iuj ri - ia t + /x 2 ) - m 2 (iu n + ^i){iuj Tl - iv m + /i 3 ) 

+m 1 (iu n -iai + fi 2 )(iu n -iu m + fi 3 ) (53) 



and for which the abbreviations E\ = y(p + Tnf, E 2 = \J(q — p) 2 + m 2 and E : 
y (q — k) 2 + m\ have been introduced. 
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As was done in the case of the polarization, it is useful to make a decomposition 
of this function in terms of elementary integrals 



T(iu m ,k;iai,p) 



N, 



57T 



(m 3 - m 2 )B (k - p, m 2 , fi 2 , m 3 , /i 3 , iv m - ia t ) 
+(m 1 - m 2 )B (p, mi, /ii, m 2 , \i 2 , ia{) 
+(777,1 + m 3 )B (k, mi, ^i, m 3 , // 3 , iv m ) 

+ m\{m 3 — 777,2) + m 2 (mi + m 3 ) + 7713(777,1 — 7712) — 2mim 2 m 3 
+m 3 (p 2 - (iai - fi 2 + /ii) 2 ) - m 2 (£ 2 - {iv m - fi 3 + /Ji) 2 
+771! ((p - A;) 2 - (ia, - iv m - fi 2 + fi 3 ) 2 

xC (p,k,mi,fii,m 2 ,fi 2 ,iai,m 3 ,fi 3 ,iu m ) . (54) 



The function Bq has already been given in Eq. (0). It is necessary to introduce a 
third loop integral Cq (cf. Appendix [AD, which is explicitly given as 



C (p, k, mi, /ii, m 2 , fj, 2 , ia h m 3 , /x 3 , iv m ) 
1 1 



16tt 2 



\P\<A 



d 3 q 
(2tt) j 



(55) 



x 



((«*;„ + /ii) 2 - £ 2 ) ((iw n - + /i 2 ) 2 - £f ) ((iw„ - zV m + /i 3 ) 2 - £f) 



As with the polarization function that was given in in Eq. (jTTj) , Eqs. (|50| 
defined intentionally without any flavor factors. 



are 



2. t and u channels 



The t and u channel exchange diagrams shown in Fig. [7| have the general form 



iM t = ft5 cl c 2 v{p2)i'lS9i) 



j>i - j> 3 - m 



(t) 



ils{igi)u{pi 



. 9192ft 
\ 

t — m' 1 ' 2 



$ Cl c 2 v (P2)75(l*i -fa + m {t) )l5u(p 



(56) 



iM u = %■ 



gmfu 

u — m^ 2 



(57) 



where, once again gi, g 2 are the quark-meson couplings for the outgoing mesons, 
and ft, f u account for flavor factors from the Gell-Mann matrices. The momenta of 
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the incoming quark and antiquark are p\ and P2, those of the outgoing mesons p% 
and P4. The mass of the exchanged fermion is denoted by and respectively. 
Note that these masses are not necessarily equal. 

For any actual calculation of cross sections and transition rates, one needs to 
sum the relevant amplitudes, take the absolute value squared and to sum over final 
and to average over initial states. Our results are listed in Appendix 0. 



C. Calculation of Cross Sections 



Having provided the prerequisites for calculation, we illustrate the main features 
via examples, (i) The process ud — > 7t + tt° provides an example in which a u channel 
diagram is required. Such a diagram occurs not only when the outgoing mesons are 
identical, (ii) The process us — ► tt + K° is chosen because it has unequal masses for 
the incoming quarks, the virtual quarks of the vertex part, and the outgoing mesons, 
(iii) The process uu — > 7r + 7r~ demonstrates the usage of mixing propagators for the 
virtual scalar mesons in the s channel. Together, these three examples indicate how 
the calculation is performed in general. 

We assume that the centre of mass system of the incoming quarks is at rest 
relative to the medium. All quantities, such as couplings g nq q, are calculated within 
this framework, and this is used in all processes calculated, e. g. in triangle diagrams. 



Since this implies that k = in Eqs. (|50|)-(l55), the vertex function T then depends 
only on the absolute value of the meson momentum. This leads to the fact that the 
total cross section depends only on the invariant energy s and the temperature T 



In what follows, we drop the momentum arguments of the scalar meson propa- 
gator T> and the three meson vertex T for simplicity. 



1. Calculation of a u ^^ +n o 

Explicit graphs for this process are shown in Fig. ^|. We have two graphs of the 
s channel type, which we label s and s', respectively. The flavor factors for these 
graphs are found to be 

fs = ~fs> = -2V2 (58a) 
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ft = -fu = -V2 



(58b) 



The relative sign of / s , f s > and f t , f u arises from the flavor matrix A3 that occurs at 
the 7T° vertex, whose uu and dd components have opposite sign. 

The virtual meson exchanged in this case is a charged a n , so that one may 
identify 

2K{ 



(59) 



1 - 4*711^, 0) 

The contribution of the vertex graph in the centre of mass system turns out in this 
case to be the same for both s channel type graphs. From our general formula 
Eq. (|51|), one has 



rv 



N r m 



c""q 



2^0(0, m q , fj, q , m q , /ig, y/s) 



(60) 



1 



+(s - 2ml)C (p 3 , 0, m q , ji q , m q , fi q , -yfs, m q , //„ yfs) 



where |p 3 | = \ Js — Ami * s ^ ne momentum of the outgoing 7r + . The total contribu- 
tion of the s channel type graphs is thus 



i{M s + M S ') = -wip 2 )uip 1 )8 ClC2 gl qq (f s V s r s + f s >V s ,T s 
= , 



(61) 



which resembles the result for the two flavor model fll2| , where this process has no 
s channel due to the lack of charged scalar resonances. 

The squared invariant amplitude thus arises from the t and u channels alone, 
and (|B6[), this takes the form 



From Eqs. 
1 



AN 2 

^ ' c s,c 



J2\M t + M u \ 



2\2 



s{m 2 q -t) — (t — ml — m q ) 

a; _ (t - miY 

s{m 2 q — u) — (u — m\ — m q ) 2 



(62) 



-2 



{u — m 2 ) 2 
(t — ml — m 2 q ){m\ + m 2 — u) + sm\ 



{t — m 2 ){u — m 2 ) 

In the limit m n 0, this reduces to an extremely simple expression. This differs 



from that given in Ref. ||12|| , in which the u channel has been omitted. 
We obtain the differential cross section for this process to be 
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da, 



ud— >7r+7T° 



1 



J2\M t + M u \ 



(63) 



dt 16?rs(s - Am 2 q ) AN* _ 

The total cross section is now constructed from the integration over the ex- 



changed momentum t as [12 



cr «J-^7r+7ro(' S ; T) 



dt 



da 







dt 



l + M^)) 2 



(64) 



The Bose distribution function fs(x) = l/(exp((3x) — 1) is introduced, since the 
presence of mesons in the heat bath leads to an enhancement of the meson creation 
process in the medium. 



2. Calculation of a uS _, n + K o 



The Feynman graphs for this process are shown in Fig. [L0[ Here only one s 
channel diagram is possible. From the graphs, we find the flavor factors 



(65a) 



ft = 2 . (65b) 
The virtual scalar meson exchanged in the s channel is a ax, so that we have 



V 



2K: 



l-4fcT 4 -nS(Vi,0) 

The contribution of the vertex part can be derived from our general formula Eq. 
to be 



(66) 



N, 



in 



m s - m q )B (p A} m q , fi q , m s , (i 8 , E 4 ) 

+ (m s + m q )B o (0, m q , /i q , m s , fi s , y/s) 
m s ml + m q (m 2 K - s + 2£ 3 (/i s - /j, q ) 

xC (p 3 , 0, m q , n q , m q , n q , E 3 , m s , /i s , ^/s) 



(67) 



where (E 3 ,p 3 ) and (i? 4 ,p 4 ) denote the four momenta of the 7r + and the K°, respec- 
tively. Because the masses of the outgoing particles are different, the momentum of 
the 7r + say, now has to be determined via 
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|p 3 | = 7T7=\/ ( s - + m K ) 2 )(s - (m n - m^) 2 ) . (68) 
From the formulae in Appendix [B], one obtains the squared transition amplitude 

as 

-L £ + M| 2 = 8 '^" f f |OT| 2 ( S - K + -,) 2 ) (69) 

'*-'*(! S.C -'"C 



K(pr) / 



■- j ( - + m 2 ) 

t — m z q \ 



+(t -m 2 - ml)(m q + m s )j 

+ 4^,^2)2 (( m - + m S " *) (* " m K + m2 s) 
+ (m 2 — t)(s — m 2 q — m 2 s ) — 2r?vn/,,/// 



with 3? denoting the real part of the function following. 

Due to the different masses of the u and s quarks, the differential cross section 
is now given as 

Again, the momentum \px\ has to be calculated using Eq. (|58]) , where m,,., have 
to be replaced by m q , m s , respectively. 
The total cross section is defined to be 

a u - s ^ +K o(s,T) = J dt dau ^ +K ° (1 + f B (E 3 ))(l + f B (E 4 )) (71) 

in this case. 



3. Calculation of a. 



UU — >7V~ r 7V 



The Feynman graphs for this process are shown in Fig. [11]. Once again, we have 
s and s' channel graphs. Since the incoming quarks have the same flavor, we have 
to consider mixing for the virtual mesons. The propagator for the s graph reads 

V s = M a „ + ~ (2M 00 + 2V2M 08 + M 88 ) , (72) 

while the propagator for the s' graph is 
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V s , = -M^ + - (2M 00 + 2V2M 08 + M 88 ) . (73) 

Here the flavor factors for the virtual meson vertices have already been taken into 
account, so that we only have to consider the flavor factors from the outgoing mesons, 
which contribute a factor \/2 at every vertex: 

f s = fs' = ft = 2 . (74) 

The meson vertex part of the s and s' graph is the same as for ud — > 7r + 7r°, as was 
given in Eq. fl60|). Since we need the sum M s + M s ', we find that the contribution 
of the ov in this combination cancels and we obtain 

- % {M s + M S ') = -iv{p 2 )u{p x )8 clC2 gl q ^ (2M 00 + 2V2M 08 + M 88 ) r . (75) 

The rest of the calculation proceeds as in the previous two cases discussed. 



IV. NUMERICAL RESULTS 



In this section, we present our numerical results. Since we are interested in the 
cross sections as an input to transport equations, we choose to give the transition 
rates w(s,T) instead of the integrated cross sections a(s,T). These two quantities 
are related by |L8| 



w{s,T) = \v Iel \a(s,T) , (76) 
where the relative velocity of the incoming particles in the centre of mass system is 



. . v (P1P2) 2 - (mxms) 2 

Kl| = - r—, . (77) 

For equal masses of the incoming quarks, m\ = rri2 = m, this can be expressed in 
terms of the centre of mass energy s as 



/ Am 2 

\v M \ = 2^1-—- . (78) 

The introduction of the transition rate as opposed to the cross sections has a fur- 
ther advantage in that it suppresses the kinematic singularity which appears at the 
threshold for the exothermic reactions qq — > MiM 2 where 2m q > rriMi + ^m 2 - 
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To commence, we choose to discuss one simple process, and decompose the calcu- 
lation into its constituent s, t and u channel graphs, in order to indicate the relative 
magnitude of these contributions. We pick the process ud — > n + r], that by Eq. ( |48b| ) 
can also be regarded as a component of the a ua calculation. The appropriate Feyn- 
man diagrams can be obtained from Fig. |^ on replacing the tt° by 77 in this figure. 
There are two s channel, one t channel and one u channel exchanges possible. In 
this case, we find f s = f s > and ft = f u , so the s channels do not cancel. We display 
these contributions explicitly in Fig. 0. Here we plot the transition rate that orig- 
inates from \M. t \ 2 (which is equal to the contribution from |.M U | 2 ), |.M S | 2 , |A^s'| 2 , 
the sum \M t + M u \ 2 , and the sum of all four graphs, \M S + M s > + M t + M u \ 2 . From 
this figure, one can see that the t and u channels give the dominant contributions. 
Although the velocity factor in Eq. (|76|) goes to zero at threshold, the t channel con- 
tribution stays finite. However, it cancels against the u channel contribution, thus 
yielding a zero total rate here. This is a rather general feature that we wish to point 
out. Such a cancellation occurs for all processes where the flavor factors of the t and 
u channel graphs are equal, f t = f u - We notice also that the s channel contributions 
are comparatively small, except in the region ^fs ~ 0.5-1 GeV, where they exhibit a 
resonance structure due to the a n exchange. At threshold, their contributions go to 
zero, since they are proportional to s — 4m 2 , as can be seen explicitly from Eq. (|BI|). 
This demonstrates the generic features that are empirically observed in all of our 
calculations: in general the t and u channels dominate, while the s channel only 
plays a role over a small region if a resonance is excited. 

We turn now to a discussion of the four possible hadronization rates, that stem 
from the four independent cross section classes cx u j, a uS , a U u and a sS that were listed 
in Table |. 



A. ud — ► hadrons 



Figure [13] shows the hadronization rate of ud — > hadrons at T = as a function 
of ^/s. Here we have shown the contributions from individual processes, in this case 
being ud — > 7r + 7r°, K + K° and ir + r]. The contribution from the process ud — > tt + 7]' 
is negligible and has thus not been shown. One recognizes from this figure that the 
transition rate is dominated by the process ud — > 7r + 7r°, as is expected physically. A 
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further physical feature that is common to all the calculations is also illustrated here, 
and we thus comment at this point: that is, that an exothermic process displays 
a finite rate at threshold. An exception to this rule follows for processes such as 
ud — > 7T + r] (which are also exothermic), where equal flavor factors occur in both the 
t and u channel exchanges and the amplitudes cancel, as was discussed above. By 
contrast, the process ud — > K + K° is endothermic, which means that the threshold 
shifts from yfs = 2m q to y/s = 2m#, where it yields a zero rate. Finally we 
comment that the relative magnitude of the three rates contributing to the total 
rate essentially reflects the available phase space. 

In Fig. |TJ], we show the hadronization rate at various values of the temperature, 
T = 0, T = 150 MeV, T = 190 MeV and T = 250 MeV as a function of y/s. The 
threshold value of ^/s at which hadronization sets in is different for the four cases 
due to the temperature dependence of the quark and meson masses. At T = and 
T = 150 MeV the dominant process, ud — ► 7r + 7r°, is exothermic, thus giving rise to 
a finite rate at the threshold y/s = 2m q . At T = 190 MeV and T = 250 MeV on 
the other hand, the pions are more massive than the quarks, so that this process 
becomes endothermic and the threshold is given in this y/s = 2m n . Thus 

the finite values at thresholds observed in the rates at lower temperatures change 
to become maxima slightly above threshold, and a zero rate is obtained at these 
thresholds. 



In Fig. [15], we show an example for the differential cross section da/d cos 9 at 
zero temperature and at y/s = 1.5 GeV. Again the contribution of ud — > 7i + r)' is 
negligible and is thus not shown. One notices a forward-backward asymmetry in 
this plot, which is solely due to kaon production. Physically one can understand 
this as follows: after the reaction ud — * K + K° has occurred, the incoming u quark 
has become a constituent of the K + , so that the produced K + preferentially takes 
the direction of the incoming u. In the other reactions, ud — > 7r + 7r°, 7i + r] and 7T + r]' 
however, the u can become a constituent of any one of the produced particles, which 
makes these reactions forward-backward symmetric. 
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B. us — * hadrons 



Figure |X6] gives the transition rate for us — > hadrons at T = as a function of 
y/s, decomposed into its constituent processes. The dominant contributions to the 
transition rate come from the processes us — > tt + K° and us — > tt°K + , which give 
the same up to a factor of two. We again note that these are exothermic processes, 
thus leading to a finite transition rate at the threshold y/s = m q + m s . For y/s 
slightly above threshold, a maximum appears in the transition rate of us — > tt + K°, 
7c°K + , that corresponds to an excitation of the ok resonance. This maximum is not 
visible in the overall summed rate at this temperature, since it is dominated rather 
by us — > K + r] production. 

Figure [L7| shows the transition rates at the temperatures T = 0, T = 150 MeV, 
T = 190 MeV, T = 250 MeV as a function of y/s. The <jk excitation, that is not 
visible at T = 0, is now clearly seen at T = 150 MeV and T = 190 MeV. At 
T = and T = 150 MeV, the threshold is determined as y/s = m q + m s , and the 
reaction is exothermic, while for T = 190 MeV and T = 250 MeV, the reactions are 
endothermic, with threshold value y/s = m n + rriK- Correspondingly, in the exother- 
mic regime, the transition rate is finite at threshold, while at higher temperatures, 
where the outgoing mesons have become more massive than the incoming quark and 
antiquark, the transition rate vanishes at the threshold value. 



C. uu — > hadrons 



The contributions of the various processes for uu — > hadrons at T = are 



shown in Fig. || as a function of y/s. Here we plot only the transition rates for the 
production of 7r + 7r~, 7r°7r°, K + K~ and tc i]. All other channels are negligible. One 
notices that the transition rate is dominated by 7r + 7r~ production. For y/s > 1 GeV, 
the K + K~ production appears to be the next most important process, followed 
by 7r°7r° and tt°i]. We observe that the contribution of uu — > 7r°7r° is rather small 
in comparision with uu — > 7r + 7r~. We can trace this back to the allowed Feynman 
diagrams and the associated flavor algebra. It turns out that the s channel exchanges 
for both processes give the same contribution in total, although they have differing 
flavor factors associated with them. On the other hand, the t channel exchange 
for 7r°7r° contributes one half that obtained for 7r + 7r~ hadronization, due to the 
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associated flavor factors (f t = 1 for uu — > 7r°7r°, f t = 2 for uu — > 7r + 7r~). In 
addition, the uu —>■ 7r°7r° also has an u channel exchange available. This tends to 
cancel the t channel exchange. Finally, the t integration runs only over half of the 
available phase space for the process uu — > 7r°7r° |]I2" |. These features result in a 



considerable reduction in the production of 7r°7r° over 7r + 7r~. We also point out that 
the production of K°K° is negligible in contrast to the production of K + K~ . This 
is due to the fact that uu —>■ K°K° proceeds only via a s channel graph, which gives 
small contributions compared to the t and u channels, as has already been discussed. 
This is a direct simulation of the Zweig rule in the NJL model. 



One final comment to the structures observed in Fig. [18] is with regard to the 
cusp seen at threshold. This is due to the excitation of the a mode with m a « 2m q . 
At t/s ~ 1.2 GeV, one notices a very weak second maximum: this comes from the 
a', appearing in the processes uu — > K + K~ and uu — > K°K°. 

Figure ^ gives the rate for uu — > hadrons at T = 0, T = 150 MeV, T = 190 MeV 
and T = 250 MeV as a function of y/s. The cusp at the threshold value of a/s at 
T = and T = 150 MeV develops to become a sharply pronounced peak. This 
is a reflection of the exothermic nature of the lower temperature processes, which 
have threshold ^/s = 2m q , as opposed to the endothermic processes, with threshold 
t/s = 2m v . In this figure, one can see the a and a' maxima and the shifting of the 
threshold due to the temperature dependence of the masses. 

D. ss — > hadrons 

Finally we show the results for ss — > hadrons for different temperatures as a 



function of y/s in Figs. 2C and 21. Here the dominant contribution comes from 
the processes ss — > K°K° and ss — > K + K~, which have identical cross sections, 
and ss — > tjt). The pion production rate is very small, since these processes do not 
have t and u channels available, again simulating their suppression by the Zweig 
rule. One notices a sharp resonance in these curves near threshold. This is due to 
the a' intermediate state. The threshold values of y/s are determined at T = and 
T = 150 MeV by \fs = 2m s , where the reaction is exothermic. At T = 190 MeV and 
T = 250 MeV, this is rather given by v^s = 2m^, since the reaction is endothermic. 
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E. Averaged Transition Rates and Hadronization Times 



In order to make a connection with transport theory, we introduce the energy 



averaged transition rate |L8 



s,(T) =MT)Wi!S 12 ^^^ • (79) 

which is the average transition rate for quarks and antiquarks coming from a thermal 
medium, their distribution being given by the Fermi function /f(x) = 1/ (exp(j3x) + 
1). The quark density is given as an integral over this function 

d 3 p 
(2^ 



Pi( T ) = I T^-^cfF ( yP 2 + m? ) (80) 



for a given quark species. In Eqs. ^), the factor 2N C accounts for the number 
of degrees of freedom. To perform the integration, we have made the approximation 
that cr(s, T) depends only on s, even when we are not in the centre of mass system. 
With this assumption, w(T) can be cast into the form 

w(T)=T dsJ(p lP2 y - ( mi m 2 ) 2 a(s,T)P(s,T) , (81) 

J (m 1 +m 2 ) 2 

where the weight function P(s,T) is given as 

P(s,T) = - 77 J-___L dEi[2N c f F (Ex)] (82) 



Pl (T)p 2 (T) 

POO 

x / dE 2 [2NJ F (E 2 )}Q(A\p 1 \ 2 \p 2 \ 2 - (s - (m\ + m|) - 2E 1 E 2 ) 

J mo 



im,2 

We show the results of the average rate w(T) for the four hadronization classes 



ud, us, uu and function of temperature in Fig. E2l One observes that the 

quantity w(T) stays fairly constant in the region from T w 50MeV up to the pion 
Mott temperature Tm-k- At Tmk, it displays a minimum. This behavior differs 
strikingly from the curves shown in Ref. [Q: whereas the cross sections shown 
there decrease towards the critical temperature and then diverge sharply, we obtain 
a minimum at the pion Mott temperature with no divergence at all. There are 
two reasons for this behaviour: (i) The reason for the occurrence of the divergence 
has been explained in Ref. |T2[ as originating from the term da/dt ~ l/(t 



m 



that arises from the t channel amplitude. These authors have worked in the chiral 
limit, so that at the critical temperature m q = 0, and the t integration required 
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for the total cross section, yields a logarithmic divergence. Here, on the other hand 
we consider the physical situation of finite current quark masses, so that the quark 
propagators stay finite for all temperatures, and no divergence may occur, (ii) The 
behaviour of the cross sections is also determined by the quark-meson couplings, 
that multiply the transition amplitude squared. In the chiral limit, as calculated in 
Ref. ||12|| , these tend to a constant value at the transition temperature, so that the 



overall behavior is given by the mechanism described in (i). On the other hand, in 
our case, when tuq ^ 0, the quark-meson couplings strictly approach zero at the 
Mott temperatures, and are thus responsible for the observed behavior. 

The little dip observed in the curves at T = 180 MeV corresponds to the Mott 
temperature of the i], Tm v - One also notes a discontinuity in the production rate for 
the class uu — > hadrons at T = 165 MeV. This occurs in the model since the a in 
the SU(3) x SU(3) case is weakly bound at T = 0, and becomes a resonance at the 



Mott temperature T Mcr = 165 MeV. This leads to a transition rate ~ 1/ — 4m^ 
and thus to a discontinuous averaged transition rate. Due to kinematic reasons, 
this is only visible in the processes uu — > 7r + 7r~ and uu —>■ 7r°7r° that contribute to 
uu — > hadrons. Note that this structure would not appear for the SU (2) x SU (2) 
NJL model, where the a is a resonant state for all temperatures. This is thus a 
model dependent feature. 

From w(T) we compute the hadronization times 

T u 1 ( T ) = w U u{T)pu + w u d(T)pj + w uS (T)p s (83) 



T s ( T ) = Wsu(T)pu + w s j(T)pj + w S s(T)p s . (84) 

The quantity t u (T) represents the mean lifetime of a u quark in a plasma before 
it hadronizes into a meson. Analogously t s (T) is the mean lifetime of a s quark. 
The numerical results are shown in Fig. [23|. In the temperature region of interest, 
150 MeV <T< 250 MeV, the hadronization times are 2-3 fm/c for a u (or d) quark, 
and 3-4 fm/c for a s quark. As expected, the hadronization time becomes infinite 
at the Mott temperature, since all hadronization cross sections go to zero here. 
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F. Strangeness Production 



In this section, we discuss the production rates for the processes that change the 
total number of s plus s quarks. The enhancement of this quantity observed in the 
Kj tt ratio in nucleus-nucleus collisions over N+N collisions is one of the observables 
which indicates new physics The processes of interest are (i) ud —>■ K + K°, 

(ii) uu — > K + K~ , K°K° and (iii) ss — > 7r + 7r~, tt ti , rjrj, rjrj', rj'rf '. While the first 
two reactions increase the total number of s + s, the third one decreases this num- 
ber. The production/reduction rates w(T) are shown in Fig. |23], respectively as a 
function of temperature. Both increasing and decreasing processes show qualita- 
tively the same behavior. The strangeness increasing processes are endothermic and 
thus occur preferentially at higher temperatures, whereas the strangeness decreasing 
processes (at least the dominant ones) are exothermic and thus can also happen at 
low temperatures. Since the processes in these three cases do not dominate their 
respective hadronization class, the overall magnitude of the rates is rather small 



compared to that shown in Fig. |22 

It is useful to define the ratio 
(AK) 

gain ^ud-^K+K^ w uu^K+K~ 4" w uu^K°K 



(85) 

measuring the total hadronization from up and down quarks into kaons relative 
to the pion production due to hadronization. Here we have assumed that the f] 
decays into three pions and we have neglected the contribution of the rf . The 
temperature dependence of this quantity is depicted in Fig. ^5|. The discontinuity 
at T = 165 MeV is due to pion production, as has been discussed previously. Kinks 
occur at the respective Mott temperatures T Ma , T Mr} , T MK and T MlT . Numerically at 
T = 150 MeV, the strangeness enhancement ratio rises to {AK) za _\ n /ix = 0.01. This 
enhancement of about 1% is to be compared with the enhancement of 9% for S + S 
collisions at 200 GeV/A over N + N collisions, as found by the NA35 collaboration 
flEfl . The loss term, defined as 

(AiT) l oss 



7T 

PsPsi^ss— *7r+7T — 4~ ^ss— >7r°7r° 4~ ^ss-^r]r}' 4~ ^ ss—*'qr) 4~ ^ ss— nj' '77') 



6) 
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is also indicated in Fig. |25|. The maximal loss is an order of magnitude less 
than the gain in strangeness due to hadronization. Note that the back reaction 
ss — > nonstrange hadrons has been calculated under the assumption of thermal 
equilibrium values for p s and p s - Should thermal equilibrium for strange quarks 
however not be reached, one would expect this loss term to be even more strongly 
suppressed. 

V. SUMMARY AND CONCLUSIONS 

In this paper, we have used the SU(3) x SU(3) NJL model to calculate tem- 
perature dependent cross sections and transition rates in a scheme ordered in the 
inverse number of colors, 1/N C . To this end, we have included an exact calculation of 
the pseudoscalar and scalar meson masses, treating these particles as resonances at 
temperatures beyond their respective Mott temperatures, when they may dissociate 
into their constituent quark and antiquark. Our calculation has been performed for 
finite values of the current quark masses, which is regarded as the physical situation. 

In describing the hadronization procedure, we have identified four classes ud, 
us, uu and ss, determined by the initial incoming quark and antiquark and final 
states 7T7T, iiK, KK, Krj and rjrj. The individual processes contributing to each class 
have been listed and calculated, and their contributions to the cross sections have 
been studied as a function of temperature and centre of mass energy. We are able 
to account for all the observed structures, that result as a function of (i) threshold 
effects, (ii) an intermediate resonance being excited, or via the (iii) summation of 
several processes. The relevance of each process is discussed, and we observe, for 
example, a forward-backward asymmetry in the class ud — > hadrons, that is due to 
the process ud — > K + K°. Averaged and total production rates are calculated for 
each class. 

What can be learned from this calculation? From the measured spectra of the 
hadrons produced in high energy heavy ion collisions, one derives temperatures in the 



range of 150-200 MeV jig] . This temperature range is also found for the location 
of the phase transition in lattice gauge calculations. In our calculation, the pion 
Mott temperature Tm-k = 212 MeV plays the role of separating the regions where 
pions are stable (T < Tm-k) and where they exist as resonances. The unphysical 
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aspect of our model is the appearance of constituent quarks for T < Tmw Thus all 
hadronization cross sections for T <C Tm* are only of academic interest. If we focus 
on the temperature region between 150 MeV and 250 MeV, two (in our opinion) 
physically relevant numbers can be extracted. 

(i) The hadronization time: The quark-gluon plasma (if it is reached at all) is 
a transient state because of the rapid expansion of the system. Therefore it is of 
great interest to know the times for the various stages of the system: thermalization 
of the quark-gluon plasma, hadronization and final state interaction time. In our 
calculation we find an average time for the hadronization of a quark of 

r had = 2-4fm/c . (87) 

This value may be considered not unrealistic in view of expansion scenarios. 

(ii) Strangeness enhancement: In experiments, a sizeable value of strangeness 
enhancement AKj tt has been observed in heavy ion collisions over N + N ones. 
This value is AK/tt = 9%. Our calculation yields an enhancement of about 1% 
as coming from the hadronization stage. This result indicates that the dominant 
contribution to strangeness enhancement must occur before or after hadronization. 

On the theoretical side, the present paper shows that hadronization cross sections 
(as other quantities too) are nonperturbative in nature and have a rather strong 
temperature dependence. Thus realistic calculations (cascade or other ones) have 
to take these aspects into consideration, if they want to be realistic. 
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APPENDIX A: LOOP INTEGRALS 

The decomposition of Feynman diagrams into elementary integrals is a well 



known technique for zero temperature [21]. At finite temperature, however, no 
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general theory has to date been given. Large SU (3) calculations are difficult with- 
out implementing a technique, that enables one to decompose all expressions into 
a few integrals only, the calculation of each of which can be done once. The eval- 
uation of these integrals for the general case is a tedious task, which will form the 
subject of a separate publication |22|. To be concise, we simply illustrate the cal- 



culation of the elementary integrals for some special cases here, confining ourselves 
to the situation of zero chemical potentials and equal masses for all fermions. This 
is already a sufficient basis for most of the calculations within the two flavor sector. 
A full calculation, including differing chemical potentials, quark masses, and arbi- 
trary kinematics (applicable in Co) can be found in a forthcoming publication ||22|| . 
Since the NJL model is nonrenormalizable, we choose a three momentum cutoff A 
as defining our regularization scheme. 

1. Computation of A 

The first loop integral has been defined in Eq. (Q) as 

, (m ,, ) ^ ? ^ / , (A1) 

|pl<A 

where E = y/p 2 + m 2 and the limit rj — > has to be taken after the Matsubara 
summation. In the case = this is easily evaluated to be 

A(m,0) = -4 / A rfp^tanh ] . (A2) 
jo E \ 2 I 

In the limit (3 — > oo, this can be evaluated analytically. For finite temperatures, the 
integral has to be performed numerically. 



2. Computation of Bq 

The second loop integral Bq is defined as the analytic continuation of 

B (k, mi, m, m 2 , fi2, iv m ) = (A3) 



]T e iuJnV J 



(3 V J (27r) 3 ((^ n + /i 1 ) 2 -E 1 2 ) ((iu n -iu m + n 2 ) 2 -EI) 

\p\<A 
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(Ei = ^p 2 + ml, E 2 = y (p — k) 2 + mf) to the real axis. The case k — 0, which we 
require for the determination of the meson masses, is singular and has to be treated 
separately. For this case, one obtains 

p 2 tanh(4f) 



£0(0, m, 0, m, 0, fc ) = 8V dp 



AE 2 - kl 



(A4) 



+iirQ(k — 2m) tanh 



1- 



2m V 



Here V denotes the Cauchy principal value of the integral. 
For finite values of \k\, B takes the form 



1 M p 
B (k,m, 0, m, 0, fc ) = -7 / dp 



tanh(f) 



log 



(P + 2p|A;|) 2 -4A; 2 £ 2 



(A; 2 - 2p\k\) 2 - Ak 2 E 2 



(AS) 



+m (&(2p\k\ - \k 2 + 2k E\) - &(2p\k\ - \k 2 - 2k E\)) 



The imaginary part of this integral can be easily evaluated analytically. Note that 
this is not covariant due to our noncovariant regularization scheme. 



3. Computation of Co 



C has been defined to be the analytic continuation of 

16tt 2 ^ „ r d 3 q 



C (p, k, mi, iii, m 2 , /i 2 , w m , m 3 , fi 3 , ia t ) = ^ e l0JnV J 



i1<A 



1 



1 



(2nY 
1 



(A6) 



x- 



((iu n + iii) 2 - E 2 ) ((iuj n - iai + fi 2 ) 2 - El) ((iu n - %v m + fi 3 ) 2 - Ei) 

In the example presented here, we not only set /i — and consider equal masses, 
but also set k — 0, \q\ = \\J s — 4m 2 , ia\ = \\fs and iv m = y/s. The result can be 
directly applied to pion production. One finds 



C o (0, q, m, 0, m, 0, v^, m, 0, \:\fs) = -, 1 =vf dp- 

^ Js(s — Am 2 ) J 



tanh (f ) 



E 



(A7) 



x 



4 /-i 

ml - C- 



yfs-2E 



log 



<7r(e(p|5i-|e + |) + e(p|9i-|e-|)) 



m 2 n - C- 



™>l -Cf 



-<7re(p|9l-|M) 
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log 



m 2 + C- 



+i7T 



tanh (^) 



s(s — 4m 2 



log 



2m! 



+ i7re(p|^I-|e + |)J 
- (s — \q\V s — 4m 2 ) 



2m 2 



's + |g| V s - 4m 2 ) 
Z7r0(|glv / s-4m 2 - |2m 2 - s|)j 



(C± = Ey/s ± £>y s — 4m 2 , £± = E^fs ± m 2 ). In the limit m, — > 0, the func- 
tions vanish identically and Eq. (|A7|) can be greatly simplified to yield the result of 



Ref. [121 for the three meson vertex. 



For a full calculation of B for arbitrary values of the quark masses and chemical 
potentials, as well as Co for the same and arbitrary kinematics, the reader is referred 
to Ref. IE2. 



APPENDIX B: SQUARED MATRIX ELEMENTS 



In this appendix, we give the technical details required for calculating the squared 
transition amplitude with unspecified general masses. One is required to average 
over initial states and and sum over final states. From the former, one obtains a 
factor l/4iV 2 . The color trace gives a factor N c . The flavor trace is accounted for 
explicitly by the factors f s , f t and /„. Our final expressions for these functions, after 
taking the spinor trace, are 

1 — . • - .» fs\9m\ 2 



AN 2 

c s,c 



2N r 



(mi + m 2 ) ) 



(Bl) 



4iV 2 



ft 1 9i92 1 



2N C (t - mW 2 ) 2 



(mj — m\ —t)(t — m 2 + m|) 



(B2) 



+ (m 



(t)2 



t)(s — ml — m 2 ) + 2m 1 m^- ) (t — m\ + m\) 



-2m2m ((t \m 



m. 



t) — 2m\m2(t + m 



(t)2> 



4iV 2 



2N C (u - m(«) 2 ) 2 

,(i»)2 



(m.4 — m i — u){u — m 2 + m 2 ) (B3) 



+(m W2 - u)(s - m 2 - mj) + 2m 1 m ( - u \u - + 

,(«)2> 



-2m2m 



(«), 



m. 



■u) — 2m\m2\u 



m v 
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for the squares of the individual channels. Here we denote the masses of the incoming 
quarks by mi, m 2 , the masses of the outgoing mesons by 777,3, m 4 an d the mass of 
of the virtual fermion in the t and u channel by my' and respectively. 
The results for the mixed terms are 

1 K A A A* fsft\gi92\ 2 £T I" 2 2 

4iVf S = 2N C t-mw4 miK ~ m4 + t} (B4) 

— r?72(?773 — 777^ — t) + mM'(s — (?77i + 777 2 ) 2 ) 



AN 2 

c s,c 



2N r 



U — 777 



(u)2 



m 2 {m\ — ml + -u) 



(B5) 



—777,! (7773 — m 2 . — m) + m^ u '(s — (ttt-x + 777,2)^) 



4iV 2 



ftf u \gi92\ 



4/V r 



x 



(f-m« 2 )(u-fflW 3 ) 



(B6) 



'mj + 777-3 — u){m\ — 777 4 + Imirn^ — u) 
+(s - (mi - 777 2 ) 2 )(m3 + 7774 - s) 

+ {m\ + 7774 — £)(mj — 777^ + 2777i777^' ) — £) 

+2m2(777i — m^)(m 2 + m 2 — t) 
+2m 2 (mi — m^){m\ + m 2 — tz) 

+2(777i — ?77 ( ^)(mi — 777 ( *))(s — (?77i + m 2 ) 2 ) 



From these formulae, it is straightforward to derive the results of Ref. ]T^[ for the 



Case 777i = 777 2 



771 



(0 



777 



W and m 3 = m 4 = 0. 
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FIGURES 



FIG. 1. The qq scattering amplitude in the random phase approximation. 

FIG. 2. Feynman diagram for the irreducible pseudoscalar polarization function. 

FIG. 3. Temperature dependence of the constituent quark masses. The solid line refers 
to the light quarks up and down, the dashed line to the strange quark. 

FIG. 4. Temperature dependence of the pseudoscalar meson masses, as well as that of 
2m q and m q + m s . Respective Mott temperatures are indicated by the solid points. 

FIG. 5. Temperature dependence of the pion (solid line) and kaon (dashed line) cou- 
pling strengths. 

FIG. 6. Temperature dependence of the scalar meson masses and 2m q . 

FIG. 7. Generic form of Feynman graphs for the hadronization amplitudes to leading 
order in 1/N C . Quarks are denoted by single lines, mesons by double ones. The three 
diagrams represent s channel, t channel and u channel exchanges, respectively. 

FIG. 8. Three meson vertex, F(iu m ,k;iai,p). 

FIG. 9. Feynman graphs for the process ud — > 7r + 7r°. 

FIG. 10. Feynman graphs for the process us — > ir + K°. 

FIG. 11. Feynman graphs for the process uu — > 7r + 7r~. 

FIG. 12. Decomposition of the rate for ud — > ir + r] at T = into its individual contri- 
butions. At the threshold, the t channel contribution is finite, the vertical line indicating 
the position of the threshold. 
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FIG. 13. Contributions of individual hadronization processes for ud — > hadrons at 
T = 0. The contribution of the process ud — > 7r + ?/ is negligible and is thus not shown. 
The rate for ud — > 7r + 7r° is finite at threshold, which is indicated by the vertical line. 

FIG. 14. Transition rates for ud — > hadrons at various temperatures. Solid line: T = 0, 
dashed line: T = 150 MeV, dotted line: T = 190 MeV, dot-dashed line: T = 250 MeV. At 
T = and T = 150 MeV, the vertical line indicates a finite rate at threshold. 

FIG. 15. Differential hadronization cross section for ud — ► hadrons at T = 0, 
v/s = 1.5 GeV. 

FIG. 16. Contributions of individual hadronization processes for us — > hadrons at 
T = 0. The contribution of the process us — > K + r/ is negligible and is thus not shown. 
The contributions of the processes us — > 7r + -ftT° and us — > 7r -RT + are finite at threshold. 

FIG. 17. Transition rates for -us — > hadrons at various temperatures. Solid line: T = 0, 
dashed line: T = 150 MeV, dotted line: T = 190 MeV, dot-dashed line: T = 250 MeV. At 
T = and T = 150 MeV, the vertical line indicates a finite rate at threshold. 

FIG. 18. Contributions of individual hadronization processes for uu — ► hadrons at 
T = 0. Only the four dominant processes are shown. The contribution of the process 
uu — > 7r + 7r _ is finite at threshold. 

FIG. 19. Transition rates for mm — ► hadrons at various temperatures. Solid line: T = 0, 
dashed line: T = 150 MeV, dotted line: T = 190 MeV, dot-dashed line: T = 250 MeV. At 
T = and T = 150 MeV, the vertical line indicates a finite rate at threshold. 

FIG. 20. Contributions of the dominant individual hadronization processes for 
ss — ► hadrons at T = 0. The contributions of the processes ss — ► K + K~ and ss — > K + K~ 
are finite at threshold, as indicated by the vertical line. 

FIG. 21. Transition rates for ss — > hadrons at various temperatures. Solid line: T = 0, 
dashed line: T = 150 MeV, dotted line: T = 190 MeV, dot-dashed line: T = 250 MeV. At 
T = and T = 150 MeV, the vertical line indicates a finite rate at threshold. 
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FIG. 22. Temperature dependence of the averaged transition rates w(T). Solid line: 
ud — > hadrons, dashed line: us — > hadrons, dotted line: ira — > hadrons, dot-dashed line: 
ss — > hadrons. 

FIG. 23. Temperature dependence of hadronization times for light quarks (solid line) 
and strange quarks (dashed line). Note the broken scale on the abscissa. 

FIG. 24. Temperature dependence of the averaged strangeness production/destruction 
rates. Solid line: ud — ► strange hadrons, dashed line: uu — ► strange hadrons, dot-dashed 
line: ss — > nonstrange hadrons. 

FIG. 25. Temperature dependence of the strangeness enhancement factor (AET)g a i n /7r 
(solid line) and the strangeness reduction factor (AK)\ oss /ir (dashed line). 
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TABLES 

TABLE I. The four independent incoming qq states and their associated outgoing two 
meson states. 



incoming 


ud 


US 


uu 


ss 


outgoing 


7T+7T 


TT+K° 


7T + 7T~ 


ir + ir~ 




K+K° 


ir°K+ 


7r°7r° 


7r°7r° 




TT + 7] 


7]K+ 


K+K- 


K+K~ 






7]'K + 


K°K° 


K°K0 








n°r] 


rjrf 








A' 


rm 








rjrj' 


rfif 








7]7] 










r/rf 
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